Adiabatic pumping through a quantum dot with coulomb interactions: 
A perturbation expansion in the tunnel coupling 
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We present a diagrammatic real-time approach to adiabatic pumping of electrons through in- 
teracting quantum dots. Performing a systematic perturbation expansion in the tunnel-coupling 
strength, we compute the charge pumped through a single-level quantum dot per pumping cycle. 
The combination of Coulomb interaction and quantum fluctuations, accounted for in contributions 
of higher order in the tunnel coupling, modifies the pumping characteristics via an interaction- 
dependent renormalization of the quantum-dot level. The latter is even responsible for the dominant 
contribution to the pumped charge when pumping via time-dependent tunnel-coupling strengths. 
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I. INTRODUCTION 

In the absence of an applied bias voltage, a mesoscopic 
conductor can sustain a DC current component if two 
or more parameters of the device (for example gate 
voltages) are periodically modulated in time. Electron 
pumping is said to be adiabatic if the parameter variation 
is slow on the scale defined by the dwell time of the 
electrons. In this case the pumped charge depends on 
the size and the shape of the pumping cycle but not 
on its detailed time evolution, i.e. it is of geometric 
naturei The idea of electron pumping dates back to 
a work of Thoulessi^ The first experiment on electron 
pumping in single electron devices was performed 
by Pothier et ali Since then much theoretical and 
experimental work has been devoted to electron pump- 
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liquid has been discussed>2^ In quantum dots pumping 
has been studied in the limit of weak interactionf2S*2& in 
the Kondo regime fSlaswell as in the Coulomb-blockade 
regimCr^ In Refs. I25l2(il the pumped charge through 
an open quantum dot is computed by means of the 
bosonization technique. Aono, in Ref. |23 uses the 
Keldysh Green's function formalism complemented by 
the assumption that the dot retarded Green's function 
takes the non interacting form (this holds true in the 
non-interacting limit as well as in the Kondo regime). 
The authors of Ref. [23 integrate numerically the master 
equation (in the Born-Markov approximation) for the 
reduced density matrix of a double-dot pump. Recently 
a general approach to pumping through interacting 
quantum dots has been put forward by relating the 
pumped charge to the instantaneous retarded Green's 
22j^2^tions of the quantum doti^Si^ 



For non-interacting systems the theory of adiabatic 
pumping, formulated by Brouwerjii is based on the 
generalization^ of the scattering approach for quan- 
tum transport to time-dependent phenomena. This 
formulation has been applied to study several aspects 
of pumping in non-interacting systems, such as the 
study of noise and decoherence^iiMi the role of discrete 
symmetries,^ the possibility of spin pumping ^^ti^si^ and 
the effect of superconducting elements and Andreev 
reflectionii^iiSiii A diagrammatic approach was used in 
Refs. 18 19 201 to calculate the pumped charge through 
a nonintcracting system using random matrix theory in 
the limit of a large number of channels in the leads. Fur- 
thermore several works investigate pumping by surface 
acoustic waves both theoretically^, and experimentally.™ 
Pumping in interacting systems has been studied much 
less so far. Quantum pumping through a Luttinger 



In this paper we study adiabatic pumping through in- 
teracting quantum dots for temperatures much above the 
Kondo temperature but much below the level spacing in 
the dot. In this case a perturbative expansion in the 
tunnel coupling between the dot and the leads is justi- 
fied. Moreover we can restrict ourselves to consider only 
one level in the dot with a strong local repulsion in the 
case of double occupancy. We aim at the understanding 
of the influence of Coulomb interaction on the pumping 
characteristics. In order to achieve this, we extend a dia- 
grammatic real-time technique'^^ that has been developed 
to describe non-equilibrium DC transport through an in- 
teracting quantum dot. As compared to the formalism in 
our recent work,^*^ the perturbative approach presented 
here, although limited to weak tunnel- coupling strengths, 
is more transparent in identifying the physical origin of 
the various contributions to the pumped charge. In par- 



ticular it is straightforward to relate the pumped current 
to the dynamics of the average charge of the dot. 

In this work we calculate the leading- and next-to- 
Icading-order contribution of the perturbation expansion 
in the tunnel coupling to the pumped charge per pump- 
ing cycle. We distinguish two cases; pumping by chang- 
ing periodically either the gate voltage and one tun- 
nel barrier or varying both left and right tunnel bar- 
rier. Considering the first case, and furthermore taking 
into account only lowest-order tunneling processes associ- 
ated with sequential tunneling, the adiabatic pump works 
analogously to a peristaltic pump. The next-order correc- 
tion turns out to be only due to a time-dependent renor- 
malization of the dot-level position that is induced by the 
combination of Coulomb interaction and tunnel coupling 
to the leads. Remarkably this effect is not masked by 
other higher-order transport processes such as cotunnel- 
ing. The situation is even more dramatic for pumping 
with the tunnel barriers. In this case, the lowest-order 
tunneling processes do not give rise to any pumping. The 
dominant pumping mechanism is, then, of higher order, 
namely pumping by making use of the time-dependent 
level renormalization. As a consequence, the gate- voltage 
dependence of the pumped charge provides a transpar- 
ent experimental access to probe quantum- fluctuation ef- 
fects. 

The paper is organized as follows: Section FlI Al intro- 
duces the model of the quantum dot. The time evolu- 
tion of the system is described by a generalized Master 
equation in section IIIBI and an adiabatic expansion is 
performed and applied to the current pumped through 
the dot in section FlI CI The expansion in the tunnel cou- 
pling is further discussed in section IIIDI The explicit 
evaluation of the formulae obtained up to here is done 
in section III El using a diagrammatic technique. In sec- 
tion ^^ results for the pumped current and the pumped 
charge are presented and discussed. 



II. MODEL AND FORMALISM 

A. Model 

We consider a single-level quantum dot with on- 
site Coulomb interaction coupled to two non-interacting 
leads. The system is described by the Hamiltonian 
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where -ff loads, -ffdot, and i^tun describe, respectively, the 
left (L) and right (R) leads, the dot, and tunneling be- 
tween dot and leads, and are given by 
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In Eqs. (^, Ccka (c^/ca) i^ t^'2 fermionic annihilation (cre- 
ation) operator for an electron with spin a —]^ j, momen- 
tum fc, energy ej, in lead a = L, R; da (dj.) is the fermionic 
annihilation (creation) operator for an electron with spin 
(7 in the dot; and n^ = d^d^ is the number operator for 
the dot electrons with spin a. The Coulomb interaction 
on the dot is described by the on-site energy U associ- 
ated with double occupation. The leads are assumed to 
be in thermal equilibrium with the same chemical poten- 
tial and to have flat bands with constant density of states 

Pa- 

By periodically changing at least two of the system 
parameters, a DC current can be pumped through the 
dot. We choose the level position of the dot e(i) and the 
tunnel matrix elements Va{t) to be time-dependent. We 
only allow for the modulus, but not the phase, of Va{t) 
to vary in time, since a time-dependent phase would cor- 
respond to a bias voltage. We define the time-dependent 
intrinsic line width r^ (i,t') = 2TTpoya{t)V*{t'), the to- 
tal intrinsic fine width T {t, t') = Tl (t, t') + Tr (t, i'), as 
well as Ta [t) = T^ {t,t) and T (t) = r(t,t). To keep 
all formulae transparent, we set h = 1 throughout the 
paper. 



B. Generalized Master equation and adiabatic 
approximation 

As described above we consider an interacting quan- 
tum dot coupled to non-interacting leads. Since the leads 
act as baths, it is convenient to trace out the degrees of 
freedom of the non-interacting lead states to arrive at an 
effective description of the reduced system. In the limit 
of temperature much lower than the level spacing of the 
dot, only one level will contribute to transport. There- 
fore the Hilbert space for the dot is four dimensional: the 
quantum dot can be empty, singly occupied with a spin- 
up or a spin-down electron or doubly occupied. These 
states, labeled by x = {0, T, i, d}, have energy Eq — 0, 

E-r = E\ 



e and i?d = 2e -t- U , respectively. In the fol- 
lowing we use a matrix notation in the four dimensional 
Hilbert space of the dot, with boldface symbols for vec- 
tors and matrices. The probabilities to find the dot in 
the respective state are p = (po,PT>-PiiPd)'^- 

The starting point of our analysis is the generalized 
Master equation for the time evolution of the system. 



J^v{t)=f dt'W{t,t')p{t') 



(2) 
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where the matrix elements W^^.x' (^' ^') '^^ ^^^ kernel 
W {t, t') describe the transition from a state x' at time 
t' to a state x at time t. For the system considered here, 
Eq. (2) defines the most general kinetic equation for 
the dot probabilities without any approximation. Off- 
diagonal matrix elements of the reduced density matrix 
for the quantum dot, that correspond to real superposi- 
tion of different states x 7^ x', do not couple to the diag- 
onal ones since x and x' differ by a conserved quantum 



number, particle number or spin. Therefore, ofF-diagonal 
elements of the reduced density matrix do not enter any 
transport quantity. Nevertheless, quantum fluctuation 
effects involving virtual intermediate states of the quan- 
tum dot in higher-order processes such as cotunneling, 
are fully taken into account by Eq. (2) by properly eval- 
uating the kernel Wi^ 

Our goal is to describe the response of the system to 
slow periodic variations of the system parameters X{t) 
with frequency O. After waiting long enough, such that 
any memory of the initial dot-state distribution p (— cxd) 
has died out, the dynamics of the system is fully de- 
termined by the explicit time dependence of the sys- 
tem parameters. The latter enters Eq. (|2} in two ways, 
namely by the kernel W {t, t') being a functional of the 
system parameters X{t) with t G [t' ,t\, and by the non- 
Markovian structure, i.e., the time derivative of p (t) at 
time t depends on p {t') at an earlier time t' at which the 
system parameters had different values. In the adiabatic 
regime it is possible to simplify considerably the form 
of the Master equation by performing an adiabatic ex- 
pansion, i.e., an expansion in the pumping frequency fl, 
assuming that fl is small as compared to both the energy 
scales that determine the decay time of the kernel and 
the time integral of the kernel (which sets the time scale 
of the system's response to the parameter's change). The 
zeroth-order, instantaneous, term of the adiabatic expan- 
sion corresponds to freezing the value of all system pa- 
rameters, X{t), at time t, which corresponds to solving 
a time-independent problem. To obtain the first-order 
correction we need to systematically collect all contribu- 
tions linear in the pumping frequency or, equivalently, 
linear in the time derivative of X at time t. For this, we 
perform a Taylor expansion of p(i') around t up to linear 
order. 
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p (i) = / dt'W [t, t') 
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The instantaneous probabilities pj are the solution of 
the time-independent problem with all parameter val- 
ues fixed at time t. They are obtained from the Master 
Eq. ((2Jl in the stationary limit. 



= W«pi^ 
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together with the normalization condition e pj = 1, 
where e — (1,1,1,1)'^, and we have introduced the 
Laplace transform 

F(z)= / dt' e- '('-'' ^ Fit -t') 



to define wf^ = wf^(z = 0+). The first adiabatic 
correction can be obtained from Eq. ||2Jl, using Eqs. Q 
and lO. We find 
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where again Wj is the Laplace transform at zero fre- 
quency and awf^ = {d/dz)w'l'\z)\^=o+. Once w['/°^ 
are evaluated and the instantaneous probabilities p^ are 
known from Eq. ®, the adiabatic corrections p^ are 
obtained from Eq. |7J together with the normalization 
condition e^pf'' = 0. 



C. Pumped charge 

The charge Q pumped in one cycle T = 1/fl is related 
to the time-dependent current Ih(t) flowing through the 
left barrier by 



Q 



T 



h{t)dt 



Furthermore, we perform an adiabatic expansion of the 
kernel W {t, t') itself. The zeroth-order term, Wj [t—t'), 
is indicated with the superscript (i) for instantaneous and 
the subscript t to emphasize that the system parameters 
X{t) — > X{t) are frozen at time t. It depends only on the 
time difference t—t' , and only parametrically on t through 
X(t). The first-order term is obtained by linearizing the 
time dependence of all parameters X{t) with respect to 
the final time t, i.e., X{t) -> X{t) -f (r - t)-^X{T)\r=t, 
and retaining only linear terms in time derivatives. This 
linear correction to the kernel is indicated by the super- 
script (a) for adiabatic. 
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Finally, we need to perform an adiabatic expansion for 
the occupation probabilities in the dot. 
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By accounting for the time evolution of the system before 
time t in a similar way as done for the dot-state prob- 
abilities above, we can express the current into the left 
lead as 



/l {t) = e I dt'e^ 

J — OO 



'W^{t,t')v{t')., 



(8) 



where W^ {t,t') = EpPW^P (i, i'), and W"P(t,i') in- 
cludes all processes associated with transitions where the 
number of electrons (with charge e) entering reservoir a 
minus the ones leaving it equals p. 

It is straightforward to perform an adiabatic expansion 
for Eq. (jHl in the same way as for the Master equation. 
The instantaneous or zeroth-order level of the adiabatic 
expansion is sufficient to describe the DC current that is 
driven through the quantum dot by an applied transport 
voltage^ It is, furthermore, sufficient for modeling rec- 
tification, i.e., the generation of a DC current component 
by applying an AC transport voltage and appropriately 



changing some system parameter in time. In the absence 
of any DC or AC transport vohage, as considered in the 
present paper, the instantaneous part of the current van- 
ishes. In order to describe pumping^ one needs to com- 
pute the first-order adiabatic correction of the current. 
Using for Eq. ((HJ the same procedure as for the Master 
equation, we find the adiabatic part of the current to be 
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D. Perturbation expansion in tunneling 



Alongside with the adiabatic expansion we perform 
a perturbation expansion in powers of the tunnel cou- 
pling strength T for both the instantaneous and the adi- 
abatic correction of the kernel W, the probabilities p, 
and the current /l (t) . We indicate the order of the per- 
turbation expansion in T by adding a superscript, i.e., 
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O(r^) for the instantaneous 
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contribution to the kernel, and similarly for Wj 
expansion of the instantaneous probabilities begins in ze- 
roth order in T, p^*^ = pf '°^ + p^''^^ + ^(r^), in order to 
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be able to fulfill the normalization condition e p. 
By expanding Eq. © in powers of F, we find that the in- 
stantaneous probabilities should fulfill the two equations, 
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together with the normalization conditions e p, 

and e'^pj =0. As discussed above, the instantaneous 
part of the current vanishes in all order in T due to the 
absence of an applied transport voltage. 

In order to determine the adiabatic corrections to the 
probabilities, we expand also Eq. |(7J) in powers of F, 
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with the normalization conditions e pj ' 
eTp(-°) 



= and 
0. We emphasize that, in order to properly 
match the powers of F on the left and right hand side 
of Eq. Hll|l , one has to start the expansion of the adia- 
batic correction of the probabilities in minus first order 
in F, pi°^ = pj"'"^^ + p^"^"^ + 0(F). At first glance, such 
an expansion might look divergent for the weak-coupling 
limit, F — > 0. However, in the validity range of the adi- 
abatic expansion everything remains well defined: the 



adiabaticity condition requires that the energy scale de- 
fined by the pumping frequency 17 is much smaller than 
the tunnel-coupling strength F. Since the time derivative 
on the left hand side of Eq. Hll|l introduces a factor fi, we 
see that p^ ' scales with fi/F, which is always much 
smaller than 1 in the adiabatic limit. 

The perturbation expansion of the adiabatically 
pumped current is derived from Eq. @ and reads 
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We see that the lowest-order contribution to the pumped 
current starts in zeroth order in the tunnel coupling 
strength F, as it consists in a product of a first- and 
a minus-first order term in the tunneling coupling, but 
scales linearly with the pumping frequency O. This con- 
trasts with the DC current driven by a finite bias voltage, 
for which the lowest-order contribution is linear in F. 

Certain properties of pumping can be derived by a 
closer inspection of the perturbative expansion of the 
Master equation. To zeroth order in F the pumped cur- 
rent is nonzero only if p(°^~i) is non vanishing, which, 
according to Eq. Hlla|l requires that the zeroth-order in- 
stantaneous probabilities p'-*'°) depends on time. How- 
ever the latter are simply determined by the Boltzmann 
factors of the corresponding state energies: 
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where E^ is the energy related to the dot state x^ 
(3 ~ 1/k^T is the inverse temperature and Z the par- 
tition function. In particular, the probabilities p^ are 
independent of the tunnel couplings. As a consequence, 
in order to have a non-vanishing zeroth-order pumped 

current /^ , one of the pumping parameters has to be 
the level position. When pumping with the two barrier 
heights, /l vanishes. 



E. Diagrammatic rules 

In order to evaluate the kernel W of the Master equa- 
tion Eq. (O we use the diagrammatic perturbation ap- 
proach to transport through interacting quantum dots 
developed in Ref. |33- While in Ref. 33 the diagram- 
matic language was derived for DC transport with time- 
independent system parameters, we generalize the ap- 
proach in this section to account for the adiabatic time 
dependence of the external parameters. 

We start with deriving the Master equation from a 
very general point of view in order to relate its kernel to 
a set of diagrams to be evaluated. In general, the (time- 
dependent) transport properties are governed by the time 
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FIG. 1: Example for the time evolution of the reduced den- 
sity matrix. The upper and lower line represent the forward 
and backward time propagation along the Keldysh contour. 
Tunneling lines connecting vertices represent tunneling events 
with the left (right) reservoir. Next to the propagators the 
respective dot states are indicated. 



evolution of the reduced density matrix of the dot ob- 
tained after tracing out the degrees of freedom of the 
non-interacting lead electrons. Since the leads are non- 
interacting and in thermal equilibrium, the lead electrons 
can be integrated out making use of Wick's theorem, i.e. 
contracting pairs of creation and annihilation operators 
'^crka (cCTfea)- Furthermore, since in the case of pumping 
there is no voltage applied between left and right lead, the 
occupation of electronic states in both leads is described 
by the same Fermi distribution function f{uj). The time 
evolution of the reduced density matrix is related to the 
propagator 11 {t, t') by 



p(i) = n(i,t')p(t') 



(13) 



Contributions to this propagator can be depicted as 
diagrams on the Keldysh contour, where contractions of 
fermion operators of the leads are indicated as tunneling 
lines. An example is shown in Fig. ^ The propagator 
n (t, t') can be expressed in terms of its irreducible part 
W (i'", t") by means of a Dyson equation: 

n (t, i') = 1 + / dt'" f dt"w (t'", t") n {t'\ t') . 

Jt' Jt' 

(14) 

The irreducible diagram part W (i"', t") is defined as the 
sum over all diagrams in which any vertical cut crosses 
at least one tunneling line (see Fig. |21 as an example) . 

Performing the time derivative of Eq. (|13|l , plugging in 
Eq. H14|l . and shifting the lower bound of the remaining 
time integral to minus infinity, we obtain the generalized 
Master equation of Eq. Q . 

In a similar way we proceed for the current /„ (t) = 
e-^{Na (t)) of particles flowing into reservoir a, which is 
given by 



/„ (i) = -iej2 [vcicU^d^m - v:{dUako.m 



(15) 
It can also be expressed in terms of diagrams. They con- 
tain a vertex at time t that has the same structure as the 
tunneling vertices in Fig. ^ and can be attached to the 
upper or the lower propagator line. We find Eq. (jS]) where 
W" (i, t') = Y.pP^"'^ (t^ t') has the following properties. 
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FIG. 2: Examples of diagrams contributing to the irreducible 
diagram part Woo- 



The terms W"^ (t, t') are given by all diagrams for which 
the number of tunneling lines with reservoir index a run- 
ning from the upper to the lower propagator minus the 
number of those with reservoir index a running in the 
opposite direction equals p. 

In the following we summarize the diagrammatic 
rules for the kernel W^ (as described in Ref. 13^ 
and discuss additional rules for the evaluation of its 
adiabatic expansion. Examples for the application 
of the rules given below are shown in appendix ^ 
We start with rules for the Laplace transform of 
Wf'"^(z) = /!^dt'exp(-z(i - t'))wf'"^ (t,t') in n-th 
order in the tunnel coupling. This will directly lead 
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US to the desired objects Wj ' = Wj ' (2:)|z=o+ and 

awf^"' = (9/9z)W^")(z)U.o+. 

(1) Draw all topologically different diagrams with n 
directed tunneling lines connecting pairs of vertices 
containing lead electron operators. Assign a reservoir 
index a, an energy to and a spin index a to each of these 
lines. Assign states x and the corresponding energies 
E^ (t) to each element of the Keldysh contour connecting 
two vertices. Furthermore, draw an external line from 
the upper leftmost beginning of a dot propagator to the 
upper rightmost end of a dot propagator that carries the 
(imaginary) energy —iz. 

(2) For each time segment between two adjacent vertices 
(independent on whether they are on the same or on 
opposite branches of the Keldysh contour) write a 
resolvent 1/AE (t) where AE (t) is the difference of 
left going minus right going energies (including energies 
of tunneling lines and the external line - the positive 
imaginary part of iz will keep all resolvents regularized). 

(3) Each vertex containing a dot operator dj gives rise 

to a matrix element {x'lda \x) where x ix') is the dot 
state entering (leaving) the vertex with respect to the 
Keldysh contour. 

(4) The contribution of a tunneling line of reservoir a is 
^Fq [t] / (cj) if the line is going backward with respect 
to the closed time path and ^Fq (t) [1 — / (w)] if it is 
going forward. 



(5) The overall prefactor is given by {—i) (—1) where 
b is the total number of vertices on the backward 
propagator and c the number of crossings of tunneling 
lines. 

(6) Integrate over the energies of tunneling lines and 
sum over reservoir and spin indices. 
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To derive the rules for the adiabatic corrections 
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we first analyze how the time-dependent pa- 



rameters enter the expression of the kernel W(i,i'). 
The time-dependent variables for which the adiabatic 
expansion has to be performed are Va (t) and e (t) , 
where the first one only appears in the product 
Tq {ti,tj) — 2TTpVa {U) Va (tj) associatcd with a tunnel- 
ing line, and the latter only in the isolated-dot propa- 
gator exp(—iJ^'dt'E^(t')] for each segment between 

adjacent vertices. While for the instantaneous kernels 
all parameters were taken at time t, now we perform a 
series expansion around the same time t and keep all 
contributions linear in a time derivative of the pumping 
parameters, 

(16) 
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The factors {ti — t) or [ti — t)^ can be included in the 
diagrammatic rules in the following way: introduce an 
additional external frequency line with the imaginary en- 
ergy —izi from the vertex at ti to the rightmost vertex at 
t (or the imaginary energy —izj from the beginning of a 
dot propagator line at tj to the rightmost upper end of a 
dot propagator line at t), performing the first derivative 
with respect to Zi (or second derivative with respect to 
Zj) then set Zi = 0+ and Zj = 0+. The external frequency 
lines are drawn as dotted lines in Fig. |3| 

The rules to compute the contribution to the adiabatic 
corrections Wf'^' due to the time-dependence of T{t) 
read: 

(7a) Add to all diagrams needed for W^ "-^(z) additional 
external frequency lines between any vertex ti and the 
right corner of the diagram and assign to them an (imag- 
inary) energy —izi. Note that an eventual external fre- 
quency line between two right corners of a diagram does 
not contribute and can always be omitted. 
(7b) Follow the rules (1) to (6) taking into account the 
extra lines. 

(7c) Perform a first derivative with respect to Zi and mul- 
tiply it by the factor i^ (i) yjj)- Sum all the contribu- 
tions obtained in this way. 
(7d) Set all the external frequencies Zi and 2; to 0+. 

The contribution to the adiabatic correction Wf"'^' 
due to the time-dependence of the level position can be 



-IZ 




-IZq 



FIG. 3: Example on how to add the external frequency lines. 
The line —iz is needed for the evaluation of 9W; —izi, —iz2, 
and —iz4 are needed for the contributions to W'^' due to 
both ^ and ^; —iz3 does not contribute and can be omit- 
ted; and —izs is additionally needed for the evaluation of the 
contribution to W'^'^' due to ^. 

at 



computed in a similar way: 

(8a) In addition to the external frequency lines added 

according to rule (7a), put one more external frequency 

line from the left corner of the diagram with no vertex 

to the right corner. 

(8b) Follow the rules (1) to (6) taking into account the 

extra lines. 

(8c) Perform a second derivative with respect to Zi and 

multiply by — ^ ^^^ ^' (t), where x ix') is the dot state 

entering (leaving) the vertex of the external frequency 
line at ti with respect to the Keldysh contour. The term 

-^ {—jf-) is omitted if the segment associated with E^ 
(Ey) does not belong to the diagram. Sum all the con- 
tributions obtained in this way. 
(8d) Set the external frequencies Zi and z to 0+. 



III. RESULTS 

In this section we show the results for the pumped 
current and the pumped charge through a single-level 
quantum dot. As no bias voltage is applied, the only 
contribution to the pumped current arises from the adi- 
abatic correction, and hence we drop the superscript (a) 
for the current. As we discuss below, the properties of 
the pumped charge, in the regime discussed in this pa- 
per, can be understood to a large extent in terms of the 
time-dependence of the occupation of the quantum dot, 
(n) = Pi + Pi + 2pd- For this reason, we first want to 
discuss the perturbation expansion of the instantaneous 
average charge occupation. The contribution to zeroth 
order in F turns out to be determined by the Boltzmann 
factors of the energies associated with the states of the 
isolated quantum dot. This yields 



(n) 



(vo) 



2/(e) 



l + f{e)~fie + U)- 



(18) 



The first-order correction accounts for quantum fluctu- 
ations due to tunneling from and to the leads. There 
are two qualitatively different effects which are due to 
tunneling and correspondingly we present the results for 
the first-order corrections as a sum of two contributions 
(n}(»>i) = („) (abroad) j^ (^)(»a-on)^ ^^lem the Contribution 
to broadening is the sum of the contributions of the two 

leads (nU^'broad) __ /^\(i, broad, R) i /^\ (i, broad, L) 

First, the dot levels acquire a finite life-time broaden- 
ing due to the coupling to lead a, which is accounted for 
by 

^^^(«,broad) ^^2- (71) ('■")) 0' (e) + (n)('^°V (e + U), 

(19) 
where 0'(w) is the derivative of (t>{uj) — 

5^Me* (i -I- ^j, and * is the digamma function. 

The first term accounts for the broadening of the 
resonance at e between empty and singly-occupied dot; 
it has a prefactor 2, when the dot is empty, and is 
zero when the dot is doubly occupied. The second 
term accounts for the broadening of the resonance at 
e + U between singly- and doubly-occupied dot (this 
contribution is zero if the dot is empty) . In the case that 
the dot is singly occupied both terms contribute with a 
prefactor 1. 

Second, the combination of tunneling and charging en- 
ergy gives rise to a renormalization of the level position, 
e — > e -I- (T (e, r, C/) with 



adiabatically pumped current. The result of the diagram- 
matic approach explained above can be written in the 
form 



a{e,T,U) = ^{e + U)-^{e) 



(20) 



as it is expected from the poor man's scaling analysis.^ 
The level renormalization is positive when the level is in 
the vicinity of the Fermi energy of the leads and negative 
for e+U being close to the Fermi energy. This means that 
the distance between the two resonances from empty to 
singly-occupied dot and from singly- to doubly-occupied 
dot is effectively decreased. The changes to the instan- 
taneous average occupation due to the level renormaliza- 
tion reads 

{nf''^'^-^^^^(^{n)^^'"^)a{e,r,U). (21) 

The sum of Eqs. (|19|l and (|21(l is directly found from 
the correction in first order F to the occupation proba- 
bility, which is shown explicitly in appendix We re- 
mark that the above interpretation of the two terms is 
in full agreement with known exact results for the non- 
interacting case. For U = (and fiat density of states in 
the leads), the level renormalization vanishes. The spec- 
tral density is then equal to the Breit-Wigner function 
and its expansion in zeroth and first order in F leads to 
the non- interacting result of (n)^*'''' and (n)(*'t>road)^ 



A. Adiabatically-pumped current 

We now proceed with solving Eqs. (|10h ,'). (jllb ') and 
(112^ ) for arbitrary interaction U to get the zeroth-order 



rf (t) = -e^#( 



\(^o) 



(22) 



This suggests the following interpretation. As the dot 
occupation is changed in time by varying the pumping 
parameters (one of them must be the level position since 
(n)*^*'°^ is independent of the tunnel-coupling strengths), 
the charge moves in and out of the quantum dot gener- 
ating a current from/into the leads. The contributions 
fiowing through barrier a split weighted by the time- 
dependent relative tunnel couplings F^/F. 

By means of Eqs. (ITUb ). ifTTb ) and ifT^ 'l. one finds the 
first-order-F contribution to the current, 



lH' it) - 



1 ^„U^,broad,L)W Iill/„)(M-en) 

dt \ J T dt 



(23) 
Again, we have written the result in such a form that an 
identification of the pumping mechanism is straightfor- 
ward. The first term of Eq. ^Z',\\ contains the contribu- 
tion due to the correction of the average dot occupation 
induced by the tunnel coupling to the left lead. Intu- 
itively the finite-life-time broadening due to the coupling 
to the left lead is associated with tunnel processes of elec- 
trons through the left barrier. Any change in the life-time 
broadening due to coupling to the left lead will, therefore, 
result in a current through the left barrier only. As a re- 
sult, this first term contains a total time derivative, and 
as parameters are periodically changing in time, it will 
not lead to a net pumped charge after the full pumping 
cycle. We conclude that changing the life-time broad- 
ening of the dot level does not contribute to adiabatic 
pumping. The second term has the same structure as the 
zeroth-order contribution, Eq. if^ . It can be understood 
as the correction term introduced by renormalizing the 
position of the dot level, which may be time dependent 
via time-dependent tunnel couplings or a time-dependent 
gate voltage. Again, the charge transferred in/out of the 
quantum dot splits into two currents to or from both 
leads with relative weight F^/F. Now, even if the dot 
level is constant and only both the tunnel couplings Fq, 
are varying in time, a finite charge can be pumped by 
means of level renormalization. 

It is useful to compare these findings with a perturba- 
tion expansion of the DC current driven by a DC trans- 
port voltage. In lowest (first) order, current is carried 
by sequential-tunneling processes. A systematic calcu- 
lation of the second-order linear conductance'^^ shows 
that quantum fiuctuation due to tunneling give rise to 
three different types of correction terms. The first one, 
which dominates the linear conductance in the Coulomb- 
blockade regime away from resonance, is due to cotunnel- 
ing. One way to depict cotunneling is to understand it as 
transport through the finite-life-time broadened dot level. 
It would, thus, correspond to the first term of Eq. I|23|l . 



Close to resonance, however, there are two more cor- 
rections to the sequential-tunneling linear conductance. 
They can be cast as sequential tunneling but with renor- 
malized level position, as discussed above, or with renor- 
nialized tunnel coupling strength. For the DC current, all 
these three contributions are present at the same time, 
which makes it challenging to identify them separately 
in an experiment. For the adiabatically pumped charge, 
where correction terms associated with a renormalization 
of the tunnel couplings and level-broadening effects van- 
ish, the situation is distinctively different. Studying adi- 
abatic pumping is, therefore, a convenient tool to access 
the energy-level renormalization. This is most dramatic 
in the case when the zeroth-order pumped current is zero, 
i.e., when pumping is done by changing both tunnel cou- 
plings. In this case, the dominant contribution to the 
pumped charge is due to time-dependent level renormal- 
ization. 




FIG. 4: Pumped charge up to first order F in units of er/i/F^ 
as a function of the time-average level position in units of F 
for different values of U. Pumping parameters are e and Fl. 
The temperature is ksT = 2T. 



B. Weak pumping 

When writing the pumped charge, we report as in- 
dices, in the following, the particular choice of pump- 
ing parameters it refers to. For example if the pumping 
fields are Fl and e, we indicate the charge as QrL,e- We 
now concentrate on weak pumping. We write the time- 
dependent parameters in the form e(i) — e + Ae(i) and 
Tait) = Ta + AT ait) (with Ae(i) and AT^it) having 
zero time average) and expand the current up to bilinear 
response in the time-dependent part of the parameters. 
Choosing e and Fl as pumping parameters we obtain up 
to first order in F: 



n ^R d 



,-Ai.o) , f, . d - (1,0) 



(24) 

where the prefactor rji characterizes the amplitudes of 
the pumping parameters as well as their relative phase: 



m 



dAe 



ATi^dt 



It was already pointed out for a noninteracting system in 
Ref. y that in the limit of weak, adiabatic pumping the 
pumped charge is proportional to the surface enclosed 
in parameter space during one pumping cycle, which in 

this case is equal to rji. Furthermore, (n) ' is the in- 
stantaneous occupation of the dot computed with the 
time-dependent parameters taken at their time-average 
value, and F = Fl -I- Fr . The first term inside the brack- 
ets is the zeroth-order-F contribution to the charge and, 
therefore, it is the dominant one. It has two peaks as 
a function of the average level position, which are lo- 
cated, in the limit U > k^T, at e = -U - kB,Th\{2) 
and e = k'QT\a.{2). The second term is first order in 
F, stems from level renormalization, and vanishes in the 
non-interacting case. The first-order-F correction tends 
to decrease the distance between the two resonances. 



In Fig. 0] the pumped charge of Eq. H24|) is shown in 
units of 6771 /F^ as a function of the time- average value 
of the level position, for different strengths of the in- 
teraction. The two peaks are directly related to tran- 
sitions between singly-occupied and empty dot and be- 
tween doubly- and singly-occupied dot. The shift by 
k-^T \n.{2) of the peak positions is of combinatoric ori- 
gin and reflects the fact that the probability of single 
occupation of the dot is increased due to the two spin 
states which lead to the same occupation number. Both 
peaks contribute with the same sign, as expected since 
the underlying pumping mechanism is the same for both 
resonances: in both cases the dot filling increases (de- 
creases) when the level position is decreased (increased). 

We now focus our attention on pumping with the two 
tunnel-coupling strengths Fl and Fr. In this case, there 
is no zeroth-order-F contribution to the pumped charge, 
as discussed above. The contribution to first order in F 
reads 



^FlTh 



^Hl (<">'■" 



■(e-,f,C/), (25) 



where 772 = /q ^ gf ^ AFRdt accounts for the pumping- 
parameter amplitudes and their relative phase as dis- 
cussed in detail before for the quantity 771 . The result for 
the pumped charge as a function of the level position is 
shown in Fig. [5| The solid line shows the result for very 
small interaction. As expected, it tends to zero, because 
the level renormalization vanishes. In the presence of in- 
teraction two peaks appear, which separate for increasing 
U . The two peaks are related to the two resonances at 
the level positions. They contribute with opposite sign. 
This reflects the opposite sign of the level renormalization 
for the two resonances. We remind that the first non- 
vanishing contribution of the perturbation expansion to 
the charge pumped through the dot by periodic change of 
the tunnel barriers is uniquely due to the effect of level 
renormalization. The height of the peaks increases for 
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FIG. 5; Pumped charge up to first order F in units of eri2/T 
as a function of the time-average level position in units of F 
for different values of U. Pumping parameters are Fl and Fr. 
The temperature is fceT = 2r. 



increasing U, growing logarithmically for large U. Even- 
tually, this increase will be cut-off by the bandwidths 
D, which we here chose to be infinite. The different 
sign of the pumped charge for the two resonances could 
serve as a signature to distinguish level-renormalization- 
induced pumping from parasitic peristaltic pumping due 
to cross capacitances of the gates modulating the tunnel 
couplings to the quantum dot. 



IV. CONCLUSIONS 

We have presented a perturbative approach in tunnel- 
ing to adiabatic pumping through interacting quantum 
dots. In particular, a general diagrammatic technique 
to perform the adiabatic expansion has been developed. 
This technique has been applied to compute the pumped 
charge through a single-level quantum dot at tempera- 
tures much higher than the Kondo temperature. Two 
pumping schemes have been considered: pumping with 
the level position and one tunnel barrier, and pumping 
with the two barriers. When pumping with the level po- 
sition and one tunnel barrier, the dominant mechanism 
of the adiabatic pump works analogously to a peristaltic 
pump. The next-order correction is related to the level 
renormalization induced by the interplay of Coulomb re- 
pulsion and electron tunneling. The situation is far more 
interesting for the case of pumping with the two barriers. 
With this pumping scheme there is no pumping in lowest 
order in the tunnel coupling, and the first non- vanishing 
contribution is due to the time-dependent level renormal- 
ization. Hence, we have demonstrated the importance of 
level-renormalization effects in pumping through inter- 
acting quantum dots. In particular, our results suggest 
that adiabatic pumping can be used to gain experimental 
access to the level renormalization in quantum dots. 
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APPENDIX A: EXAMPLES OF DIAGRAMS 

In this section we show how to apply the diagrammatic 
rules for the matrix element (Wt)o.o- 
Wc start with the instantaneous term to lowest order in 
the tunnel coupling, (W^ (z))o,o- The corresponding 
diagrams are shown in Fig. Two topologically dif- 
ferent diagrams contribute, and each of them has to be 
summed over the spin index a and the lead index a, and 
to be integrated over uj. We obtain: 



W['^'\z) 



fT,a 



duj 
2^ 



Tafiuj) Taf{uj) 



By letting z = 0+ and making use of l/(x + iO+) = 
P/x—iTTS{x), where P indicates Cauchy's principal value, 
we get 






0,0 



0,0 



-2r/ (e) , 

2r^ f ^^IM. 



n de 



U! - 



Eor the adiabatic correction we need to introduce ad- 
ditional external frequency lines according to the rules 
7a and 8a (see EigEJ. The additional line of rule 8a, 
going from the left corner of the diagram with no ver- 
tex to the right corner, does not contribute in this case 
and we have omitted drawing it. The evaluation of these 
diagrams leads us to the result 



(w(-^') 



0,0 



dt " 

TT de 



duj 

p UJ 



/H 



^ dt "■ 

■K de^ 



dio 

p LO 



/M 



We now calculate the second-order-E contribution to 
the same matrix element of the instantaneous kernel. All 
diagrams contributing to the matrix element (W^ )o,o 
are depicted in Eig.0 We sum over all appearing indices 



-IZ 



-IZ 




FIG. 6: Diagrams contributing to the first-order-F part of 
the adiabatic correction to the matrix element (Wt)^^. All 
appearing reservoir and spin indices a, a are to be summed 
over and the energy u> is to be integrated over. 
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FIG. 7: Diagrams contributing to the second-order-F instantaneous part of the matrix element (Wf)(j q. All appearing reservoir 
and spin indices a, a' , a, a' are to be summed over and the energy u is to be integrated over. 



a, a', a, a' . The variables a denotes the opposite spin of 
a. As an example, we report the result for the sum of 
the first diagram of the first and the second line (after 
setting z = 0+): 



TT 



de p uj — e 



de Jp Lu — e 



(',2) 



To obtain the full second-order contribution (Wj ' '^)o,o 
we need to evaluate the remaining diagrams in Fig. [7| 
along the same lines as discussed in this Appendix. 



APPENDIX B: OCCUPATION PROBABILITIES 

Some intermediate results were not presented in the 
main part of this article as they were lengthy or not im- 
mediately necessary for the interpretation of the pumped 
current. Here we discuss in detail the first corrections 
to the occupation probabilities. The corrections to the 
occupation probabilities are used for the evaluation of 
Eq. H12|l but do not appear directly in the results for 

the pumped current. We find for p^ ' , the first order 
adiabatic correction in minus first order in F: 



Pi 



(a,-l) 



dp. 



iU)) 



1 



1 



dt 2F [1 + / (e) - / (e + U)] ' 



The adiabatic correction of the occupation probability is 
proportional to the time-derivative of Pj , and is there- 



(i 1) 

fore an eigenvector of the matrix W^ ' . The sign of this 
correction depends on the sign of the time derivative of 
eit). 

The first-order-F correction to the instantaneous occu- 
pation probability is: 



A^,i) 



dp; 



{^,0) 



de 



-a(e,F,C/) 



1 1 
1'2'2'0 



+ (2-(n)(»'O))0'(.)(^ 



It consists of a part due to level renormalization (first 
row) and a part due to level broadening (second and third 
row). The correction due to level renormalization affects 
all four probabilities in the same functional way. 
The correction due to broadening has two contributions. 
The first one is related to the broadening due to fluctua- 
tions between empty and singly-occupied dot. It is zero in 
the case that the dot is doubly occupied and largest when 
the dot is empty. The second contribution is related to 
the broadening due to fluctuations between singly- and 
doubly-occupied dot. It is zero when the dot is empty 
and largest when the dot is doubly occupied. 
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